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Abstract Pulsar Wind Nebulae (PWNe) are bubbles or relativistic plasma that form 

when the pulsar wind is confined by the SNR or the ISM. Recent observations have 

M-J ' shown a richness of emission features that has driven a renewed interest in the the- 

P^ , oretical modehng of these objects. In recent years a MHD paradigm has been de- 

(-H ' veloped, capable of reproducing almost all of the observed properties of PWNe, 

Oj, shedding new light on many old issues. Given that PWNe are perhaps the nearest 

systems where processes related to relativistic dynamics can be investigated with 
high accuracy, a reliable model of their behavior is paramount for a correct under- 
^ ' standing of high energy astrophysics in general. I will review the present status of 

C^ , MHD models: what are the key ingredients, their successes, and open questions that 

still need further investigation. 
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When the ultra-relativistic wind from a pulsar interacts with the ambient medium, 
either the SNR or the ISM, a bubble of non-thermal relativistic particles and mag- 
netic field, known as Pulsar Wind Nebula or "Plerion" (PWN), is formed. The Crab 
Nebula is undoubtedly the best example of a PWN, and it is often considered the 
prototype of this entire class of objects, to the point that models of PWNe are, to a 
large extent, based on what is known in this single case. The first theoretical model 
^^ ' of the structure and the dynamical properties of PWNe was presented by Rees & 

H I Gun 188], further developed in more details by Kennel & Coroniti IMIISI (KC84 

hereafter), and is based on a relativistic MHD description. 
The MHD paradigm is based on three key assumptions: 
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that the larmor radii of the particles is much smaller than the typical size of the 
nebula, and particles are simply advected with the magnetic field. This is true 
up to energies of order of the pulsar's voltage, where the larmor radius becomes 
comparable with the typical size of the system. 

That radiative losses are negligible, or at least that they can be accounted for 
by renormalizing the pulsar spin-down luminosity. This again can be proved to 
be true in the case of Crab Nebula (and to some extent also in other systems 
with good spectral coverage), where the synchrotron spectrum shows that the 
particles carrying the bulk of the energy have a typical lifetime for synchrotron 
cooling longer than the age of the nebula. 

That we are dealing with almost pure pair plasma, and dispersive or hybrid effects 
(separation of scales) due to the presence of heavier ions are absent. While there 
is no direct evidence for the absence of ions, standard pulsar wind theory, and 
the success of the MHD model of PWNe suggest that, from a purely dynamical 
point of view, there is no need for this extra component. 



In it simplest form II88I the MHD model of PWNe can be summarized as follow 
(see Fig.|3]i: the ultra-relativistic pulsar wind is confined inside the slowly expanding 
SNR, and slowed down to non relativistic speeds in a strong termination shock (TS). 
At the shock the plasma is heated, the toroidal magnetic field of the wind is com- 
pressed, and particles are accelerated to high energies. These high energy particles 
and magnetic field produce a post-shock flow which expands at a non relativistic 
speed toward the edge of the nebula. 

Despite its simplicity the MHD model can explain many of the observed prop- 
erties of PWNe, and until now no observation has been presented that could rule 
it out. The presence of an under-luminous region, centered on the location of the 
pulsar, is interpreted as due to the ultra-relativistic unshocked wind. Polarization 
measures Illl2lll061l93ll57l[82ll42ll70ll53l show that emission is highly polarized 
and the nebular magnetic field is mostly toroidal, as one would expect from the com- 
pression of the pulsar wind, and it is consistent with the inferred symmetry axis of 
the system. The pressure anisotropy associated to the compressed nebular toroidal 
magnetic field ifTTlllOll . explains the elongated axisymmetric shape of many PWNe 
(i.e. Crab Nebula, 3C58). The MHD flow from the TS to the edge of the nebula also 
leads to the prediction that PWNe should appear bigger at smaller frequencies: high 
energy X-rays emitting particles are present only in the vicinity of the TS, having 
a shorter lifetime for synchrotron losses, compared to Radio particles which fill the 
entire volume, having negligible losses on the age of the nebula. This increase in 
size at smaller frequencies is observed in the Crab Nebula lll051[T4l l7ll. 

However one must bear to mind that not all properties of PWNe can be explained 
within the MHD framework, which, ultimately, only provides a description of the 
flow dynamics. For example, the acceleration of particles at the TS that accounts 
for the continuous, non-thermal, very broad-band spectrum, extending from Radio 
to X-rays ifTOSl lSl fTTOlfTTTIISSll . is usually assumed as given. The MHD model pro- 
vides no hint to the reason why the injection spectrum looks like a broken powelaw, 
with no sign of a Maxwellian component at lower energies. Moreover the MHD 
description might prove faulty if applied to particles responsible for the emission in 
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the 10-100 MeV band, whose larmor radii are comparable to the size of the TS, and 
can lead to wrong conclusion on their expected behavior. 

The MHD model of PWNe has been used, by comparing observations with the 
predictions of numerical simulations, to constrain some of the properties of the pul- 
sar wind, at least at the distance of the TS. While it is not possible to derive the 
Lorentz factor of the wind, or its multiplicity, it is possible to constrain the ratio 
between Poynting flux and kinetic energy, the latitudinal dependence of the energy 
flux, and the presence of a dissipated equatorial current sheet. This shows that neb- 
ular properties can be used to derive informations on the conditions of the pulsar 
wind at large distances. 



2 Jet-Torus structure and Inner flow properties 

Let start our discussion from young objects. The Crab Nebula, 3C58, MSH 15-52, 
G21.5 all belongs to this group. Younger objects are the most well studied and per- 
haps the ones for which the MHD models have provided the greatest insight. These 
systems are characterized by a simple interaction with the confining SNR, they are 
bright, we have broad band data, and the pulsar proper motion can be neglected. 
Older systems are often subject to a much more complex interaction with the SNR, 
they are affected by the pulsar proper motion, and usually lack the deep observa- 
tional coverage of the younger counterparts. For this reasons, models of old objects 
have also progressed far less than for young ones, and the agreement with obser- 
vations is mostly qualitative. We will leave a description of MHD models of the 
evolution of PWNe to Sec. |3] 

The KC84 model has been for a long time the reference for the understanding 
of young PWNe, with only minor theoretical developments. Things have changed 
recently thanks to high resolution optical and X-rays images from HST, CHANDRA 
and XMM-Newton, that have have shown that the properties of the emission at high 
energy cannot be explained within a simplified one dimensional model. This refers 
not just to the geometrical features that are observed, but in practice to all aspect of 
X-ray emission. 

These new data show that the inner region of young PWNe is characterized by 
a complex axisymmetric structure, generally referred as jet-torus structure (Fig.|2]l. 
First observed in Crab 154111101 . it has subsequently been detected in many other 
PWNe l50l|44l|52l[87l|45l|7l|9T]|95l[35l|9l, tothe point that the common con- 
sensus is that, with deep enough observations, it should always be detected. This 
structure is characterized by an emission torus, in what is thought to be the equa- 
torial plane of the pulsar rotation, and, possibly, a series of multiple arcs or rings, 
together with a central knot, almost coincident with the pulsar position, and one or 
two opposite jets along the polar axis, which seem to originate close to the pulsar 
itself. Even if the existence of a main torus could be qualitatively explain as a con- 
sequence of a higher equatorial energy injection 12011 it is not possible to reproduce 
quantitatively the observed luminosity. Shibata et al. l94| were the first to point that. 
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the difference in brightness between the front and back sides of the torus in Crab 
Nebula, requires a post-shock flow velocity ^ 0.4 — 0.5c, much higher than what 
expected for subsonic expanding flows. The same conclusion applies to all the other 
systems where a torus is observed. The existence of an inner ring, detached from 
the torus, and of the knot which seemed to be located inside the wind region, are 
incompatible with the assumption of a smooth flow from the TS. From a theoretical 
point of view however, the most interesting feature is the jet iTSl . because theoret- 
ical lfT2l [T3I and numerical 1341 [191 JST] |68] [30l studies of relativistic winds from 
pulsars have shown no presence of collimated energetic outflow. To this, one must 
add other observed properties, like the X-rays photon index maps of Crab Nebula, 
Vela and Kes 75 ISSllMllSSl . which hardens moving from the inner ring toward 
the main torus, while steepening due to synchrotron losses is expected, and the rel- 
atively large size of the X-ray nebula in Crab compared to the radio ||4l. The fact 
that the symmetry axis of the jet-torus corresponds to the major axis of the nebula, 
leads immediately to the conclusion that the toroidal magnetic field is paramount in 
shaping the inner flow. 




Fig. 1 From left to right: structure of the post shock flow in PWNe J40J, the funneUng of the wind 
(A) into an equatorial flow (B/C) is clearly evident. Same structure showing the complex flow 
dynamics that develops downstream of the shock and the corrugation of the shock surface which 
manifests itself as time variability at high energy (33|. Numerical result of the internal dynamics 
in the body of the PWN J40J where the flow is diverted back toward the axis by the magnetic hoop 
stresses, and is collimated into a jet. 



The keys in understanding the jet-torus structure are the magnetization and en- 
ergy distribution in the pulsar wind. It has been known for a long time lISTl . and 
has been recently confirmed with numerical simulations [,19. i68.i30l, that far from 
the Light Cylinder a higher equatorial energy flux is expected. It is this particular 
latitudinal distribution of the pulsar spin-down luminosity which naturally produces 
an oblate TS with a cusp in the polar region 12011211 . giving rise to a complex post- 
shock dynamics. The obliquity of the TS at higher latitudes, forces the flow in the 
nebula toward the equator with speeds ~ 0.3 — 0.5c. Hoop-stresses are more efficient 
in the mildly relativistic flow, and the collimation of a jet occur in the post shock 
region 117511631 . The evident complexity of this scenario makes clear that the only 
possible way to proceed require the use of efficient and robust numerical schemes 
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for relativistic MHD 1671 [39l |48l . Thanks to numerical simulations this quaUtative 
picture has been developed into a quantitative model which has been successfully 
validated against observations. 

The starting point of the MHD model is the structure of the force-free pulsar 
wind: the energy flux in the wind has a strong latitudinal dependence of the form 
L( )= Lo ( 1 + o: sin ( )) , where a is a measure of the pole-equator anisotropy, while 
the magnetic field in the wind B(0) <=<: sin(0). Various numerical simulations of the 
interaction of such wind with the SNR ejecta have been presented 1691 l40l l22l 1411 
|1091[33ll : the result of the anisotropic energy distribution in the wind is that almost 
all of the downstream plasma is deflected toward the equatorial plane, and flow 
channels with velocity ~ 0.5c can form, the value expected in order to justify the 
luminosity distribution in the torus of the Crab Nebula l94l . Shear and instabilities 
tend to destroy this collimated equatorial flow before it reaches the edge of the 
nebula, however a bulk equatorial motion survives to distances corresponding to the 
location of the torus. It is this flow inside the nebula, with speed in excess of c/3 that 
advect freshly injected particles to larger distances, giving rise to a more extended 
X-ray nebula, that the simple ID model would predictQ. The post shock flow is 
independent on the specific values of Lorentz factor or density distribution, but it 
is only determined by the puslar spin-down energy distribution, and the nebular 
dynamics cannot be used to constrain the value of the wind Lorentz factor or the 
multiplicity in the wind. 

As the flow expands away from the TS, toward the edge of the nebula, the mag- 
netization increases, until equipartition is reached. Due to the magnetic field distri- 
bution in the wind, equipartition is first reached close to the equator than at higher 
latitudes. The magnetic pressure prevents further compression beyond equipartition, 
hoop stresses in the mildly relativistic postshock flow become efficient, and the flow 
is diverted back toward the axis. This is the process that causes the formation of 
a collimated jet along the axis itself. The wind magnetization regulates the forma- 
tion and properties of the jet: for low values a < 0.001, equipartition is not reached 
inside the nebula, and no jet is formed. At higher magnetizations equipartition is 
reached in the close vicinity of the TS, and most of the plasma ends in a jet. The 
plasma speed in the jet is ^ 0.7c, in agreement with observation of the jet in Crab 
Nebula and MSH 15-52 lTT0ll56ll72l[38ll . for magnetization values a - 0.1. Asso- 
ciated with this collimated back-flow there is a global circulation inside the nebula, 
with typical speeds ^ 0. Ic that might lead to mixing with cold ions |80|. 

Numerical models offer the possibility to investigate different distributions of 
magnetic field in the wind. In particular, for oblique rotators, while the energy dis- 
tribution in the wind is identical to the aligned case llT8ll96l , the magnetic field is 
supposed to give rise to a striped equatorial region, with alternating polarities. If 
this striped wind region is dissipated, and this can happen either in the wind l79ll66l 
or at the termination shock itself 1176117711 . then a low magnetization equatorial flow 
is expected, which adds complexity to the flow structure inside the nebula. One of 
the main success of the MHD model is that emission maps based on the results 
of numerical simulations give different observational signatures if an unmagnetized 
equatorial sector, corresponding to the striped wind region, is present or not: a large 
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Fig. 2 Upper panel: left - simulated X-rays synchrotron map based on numerical MHD simulations 
of the flow |33| ; right - CHANDRA image of the central region of Crab Nebula. The two images 
have similar scaling. Lower panel: left - simulated optical synchrotron map based on numerical 
MHD simulations of the flow (different from the above one); right - HST image of the wisp region 
in Crab Nebula |56| . Note the agreement between the observed features and the results of MHD 
models. 



Striped region is needed to explain the observes inner-ring outer-torus structure of 
many PWNe, while models without it lead to single ring nebulae. 

When comparing observations with emission maps based on the fluid structure 
derived from relativistic MHD simulations we clearly see that, within the MHD 
regime, it is possible to recover almost all of the observed features, with correct size 
and luminosity. However, until now, little work has been devoted to investigate if it 
is possible to discriminate among various particle injection mechanism. Work has 
mostly focused on X-ray, and a uniform injection in the form of a single power- 
law distribution has been assumed Ii41l i109, 69i i33L Moreover one should be aware 



MHD models of Pulsar Wind Nebulae 



that even today, with better computational faciUties, it is not possible to conduct an 
exaustive sampling of the entire parameter space characterizing the interaction of 
PWNe with SNRs. Works has mostly focused on reproducing Crab Nebula, where 
many parameters are constrained by a rich set of observations. 

In Fig. |2] CHANDRA and HST images of the Crab Nebula are compared to 
maps based on a simulations with striped wind in X-rays and optical. The knot and 
the inner ring are both present, and they are due to the high velocity flow in the 
immediate post shock region, at intermediate latitudes. The main torus is visible at 
larger distances, as well as features like the anvil which corresponds to the back- 
ward side of the nebula. X-ray maps of the spectral index based on simulations also 
agree with the main observed properties of the Crab Nebula |83| and with recent 
results about Vela ||6T| and Kes 75 ||85l : in particular the spectrum appears to flatten 
moving away from the pulsar toward the main torus, without the need to assume 
any re-acceleration. All this rich emission pattern is ultimately related to Doppler 
boosting effects: at high speed, the emission is enhanced (rings) and the spectrum 
is harder. There are still problems to recover the correct luminosity /spectrum in the 
jet. This might be indicative of some form of dissipation and re-energization along 
the axis; possibly associated with local instabilities in the toroidal magnetic field 
191 which present axisymmetric simulations cannot address, but for which there are 
many observational evidences lISTl [84l l72l l38l . 

Preliminary results 1331 show that, to realistically respoduce the X-ray emission 
from Crab Nebula, a ~ 0. 1 is required (about two orders of magnitude higher than 
the ID estimate by KC84) in conjunction with a large (^ 45°) striped zone. 




Fig. 3 Variability in the wisps region in MHD 1331 . Left figure: variability at a selected point on 
axis, note the typical 2-years long cycle and also the shorter timescale variations. Right figure: 
variability of the nebula obtained by subtracting two images at 2-year distance. The outgoing wave 
pattern is easily seen. Compare with optical images from 1561 . 
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Perhaps the most promising clues to investigate the flow dynamics inside PWNe 
in the future might come from polarization, and possibly X-ray polarization: while 
emission maps mostly trace the flow velocity inside the nebula, they have little sen- 
sitivity to the presence of a small scale disordered component of the magnetic field. 
There are several indications from optical polarimetry in the Crab Nebula 1531 . as 
well as indication from recent MHD simulations lf33l 11091 . suggesting that turbu- 
lence might be present, in the body of the nebula, leading to a partial randomiza- 
tion of the field. The effect of the flow velocity on the polarization angle has been 
discussed by ||29l and BTl . and the results generally agree with available optical 
polarization measured in Crab ifTTH [T06l l93l l57l l53l . 



2.1 Time variability 

It is known that, close to the supposed location of the termination shock, PWNe 
show a short time variability mainly detected in optical and X-ray bands. This how- 
ever does not affect the global properties of the main observed features belonging 
to the jet-torus structure (i.e. the inner ring, torus and jet) which appear to be quite 
persistent on long time-scales. Variability of the wisps in the Crab Nebula has been 
known for a long time ll56l [T5l. Recent observation have shown that the jet in Vela 
appears to be strongly variables ll86l[87l . together with the main rings ||6T1 . Vari- 
ability is also observed in the jet of Crab lf84ll72l . and have recently been detected 

in MSH 15-52 IMI- 

In the strongly toroidal field of these nebulae, the jet variability, which usually 
has time-scale of years, is likely due to kink or sausage mode, or even to fire-hose 
instability 11001 . On the other hand the wisps show variability on shorter time-scales 
of months: the variability takes the form of an outgoing wave pattern, with a possible 
year-long duty cycle. 

For a long time the only model capable of reproducing the observed variability 
was the one proposed by Spitkovsky & Arons ||97ll . based on the assumption that 
ions are present in the wind. The idea of ions was also supported by kinetic simu- 
lations of acceleration in a strong shock 13]|. The presence of particles with larmor 
radii, of order of the size of termination shock, introduces kinetic effects related to 
the separation of scales, leading to compression of electrons. The model however 
requires a large fraction of pulsar spin-down energy in the ion component, which 
contrast the basic idea of leptonic dominated systems, and in general is not sup- 
ported by spectral model of PWNe. 

The most recent achievement of the MHD nebular models have been the ability 
to reproduce the observed variability ll33]|3n[Tolll09l . It is the fundamental mul- 
tidimensional nature of the problem that allows for variability of the flow pattern. 
It was already noted in early simulations 1221 that the synchrotron emissivity in- 
side the nebula varies. There were also evidences 1691 suggesting that the nebular 
flow might have a feedback action on the TS, causing it to change shape, and thus 
inducing a change in the appearance of the wisps. 
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This picture has been confirmed recently |[33lll09| : within the MHD regime it is 
possible to recover the variability, the outgoing wave patter, its typical speed, and 
luminosity variations. It is found indeed that a SASI like instability is present. Waves 
injected at the termination shock can propagate toward the axis, feeding back on the 
termination shock and triggering the injection of new waves. Simulations show that 
there is a typical duty-cycle of about 1-2 year; more generally the duty-cycle will be 
of the order of the radius of the termination shock divided by the typical propagation 
speed ^ 0.5c. 



2.2 Gamma rays 

At the moment the most promising observational avenue in PWNe research is the 
study of gamma-ray emission. HESS has shown that many extended gamma-ray 
sources are associated with PWNe both young ones in the free expansion phase and 
older ones undergoing reverberation ||36ll47l . New results from the FERMI satellite 
are just arriving HHJ], with data extending from the high energy MeV synchrotron 
part of the spectrum, to the 100 GeV IC part. 

The emission at MeV energy observed from young object should enable us to 
put constraints on the acceleration mechanism. It seems that an exponential cutoff 
in energy can explain the observed data. The more interesting question regard how 
the observed variability of the wisps manifest itself at higher energies. Naively one 
might expect that the emission at high energies should show similar variability, on 
comparable time-scales, but with amplitudes of order unity. On the other hand pre- 
liminary results 1731 suggest that in the lOOMeV range no variability is detected. 
This has important implications in term of acceleration: particles responsible for 
the lOOMeV emission have typical larmor radii of order of the size of the TS. Co- 
incidentally this is also the coherent length of the turbulence that is at the base of 
the MHD variability 1331 . In this sense high energy particles, are decoupled form 
the MHD flow, and their response to the MHD turbulence is incoherent. This of 
course deserves further investigation. In particular one would like to know at what 
energy variability reaches a maximum, where the particles start to decouple form 
the MHD flow, and to what degree this incoherent interaction with turbulence can 
induce variability. 

The emission at GeV energies is assumed to be from Inverse Compton scattering 
on background radiation, and in the case of Crab the self synchrotron. Interestingly 
the spectral properties of the comptonized radiation can be used to derive infor- 
mation about particles which are supposed to emit synchrotron in the UV, and are 
usually not directly accessible given the high UV absorption in the ISM. In the 
case of MSH 15-52 [2| Fermi results, have ruled out previous EGRET measures, 
and shown that the undetected high energy part of the particle distribution function 
must be harder than previously assumed. Similar results could be expected for Kes 
75. GeV emission could in principle help provide an independent constraint on the 
maanetization in the nebula. 
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An alternative contribution to the gamma ray emission, if protons are present 
in the pulsar wind, is the p-p scattering and related pion decay |l5]|59l. However 
uncertainties on the target number density and the spin-down energy in protons, 
make this channel hard to constrain. At present IC scattering can fit most of the 
objects, if one allows for fluctuations of the local background of order unity with 
respect to the galactic average. 

Recently a series of simplified evolutionary models for the high energy emission 
from PWNe have been developed extending beyond the free expansion ll37l l49l 
[32I . This is quite important given that the majority of the gamma-ray PWNe are 
supposed to be post-reverberation objects, where the interaction with the SNR shell 
can play a major role 

For example in many of these objects the bulk of the gamma emission is not 
centered on the pulsar, and the displacement is too large to be explained in term 
of a moving pulsar leaving behind a relic PWN. Two possible explanations have 
been invoked: an off-center compression by the reverse shock fTE\ : or the formation 
of a bow-shock tail ll28l l62l . where particle responsible for the inverse Compton 
emission in the gamma-ray band can be advected to large distances from the pulsar. 



3 Evolution of PWNe 

In the previous section we have devoted our attention to the successes of the 
fluid/MHD model regarding the emission properties observed in young system. Nu- 
merical simulations offer also a way to follow the evolution of a PWN and its inter- 
action with the SNR at later ages. 

At the moment, however, the study of old object has been limited to a qualitative 
analisys of the interaction in an attempt to recover the main phases of the evolution 
and understand how the observed multi-wavelength morphology depends on the 
interaction itself. 

In the analytic model developed by KC84 the SNR has only a passive role, pro- 
viding the confinement of the PWN. Given the complexity of the PWN-SNR interac- 
tion, a detailed study of the evolution of the system, has been possible only recently, 
thanks to the improvement in computational resources II 1021 [1611251 II 041 l43l . 

By comparing the energy in the SNR (~ 10^' ergs) to the total energy injected 
by the pulsar during its lifetime (^ 10'*'^ ergs) it is easy to realize that a PWN cannot 
significantly affect the SNR, while the evolution of the SNR can have important 
consequences for the PWN. 

There are three main phases (for a more complete discussion of PWN-SNR evo- 
lution see ||89l and gS)), in the PWN-SNR evolution. At the beginning the PWN 
expands inside the cold SN ejecta. The SN ejecta are in free expansion, so this phase 
is generally called free expansion phase. This phase lasts for about 1000-3000 yr, 
and during this period the pulsar luminosity is high and almost constant. This is 
die present phase of the Crab Nebula, 3C58, MSH 15-52, G21.5, and PWNe in this 
phase are expected to shine in high energy X-rays emission. The expansion velocity 
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Evolution of the PWN r.ontoct disooniinuitv 




Fig. 4 Left picture: schematic representation of the global structure of the PWN in the first phase of 
its evolution inside a SNR. From the center the various regions are: 1- the relativistic pulsar wind, 
2- the hot magnetized bubble responsible for the non thermal emission, 3- the free expanding ejecta 
of the SNR, 4- the ISM. These regions are separated by discontinuities: a- the wind termination 
shock, b- the contact discontinuity between the hot shocked pulsar material and the swept-up SNR 
ejecta, c-the front shock of the thin shell expanding into the ejecta, d- the reverse shock of the SNR, 
e- the contact discontinuity separating the ejecta material from the compressed ISM, f- the forward 
SNR shock. Right picture: evolution of the PWN size, from free-expansion to sedov phase (from 
(251). 



of PWNe in this early stage is typically few thousands kilometer per second. For 
this reason one can neglect the pulsar kick (velocities in the range 50-300 km/s) 
in modeling young objects, and assume the pulsar to be centrally located. As the 
system expands inside the high density, cold, supersonic ejecta of the SNR, a thin 
shell of swept-up material is formed. Given that the density of the shell is much 
higher than the enthalpy of the relativistic plasma, the shell is subject to Rayleigh- 
Taylor instability. This is supposed to be at the origin of the filamentary network 
of the Crab Nebula ||55]|60l|271, and 3C58 |17J. In the thin-shell approximation, it 
has been shown that it is possible to derive an analytic self-similar solution 1271 . 
describing the expansion of the nebula in this phase. 

The PWN will expand until eventually it will come into contact with the reverse 
shock in the SNR shell. In the absence of a central source of energy, like a pul- 
sar, the reverse shock is supposed to recede to the center of the SNR in a time of 
order of 5000-10000 yr Ii991 . From this moment on the evolution of the PWN is 
modified by the more massive and energetic SNR shell: the PWN undergoes a com- 
pression phase generally referred as reverberation phase, that can last 5000-10000 
yr. In the simple ID scenario the pressure in the compressed PWN will rise to bal- 
ance the compression and to push back the ejecta, and the nebula might undergo 
several compression and rarefaction cycles. This is however an artifact of the ID 
geometry, and conclusion based on the existence of these oscillations are not re- 
liable. More appropriate multidimensional studies have shown that the SNR-PWN 
interface is highly Rayleigh-Taylor unstable during compression |fT6l, which can 
cause efficient mixing of the pulsar wind material with the SNR. This mixing will 
most likely prevent any oscillation and the system might rapidly relax to pressure 
equilibrium. This reverberation phase is supposed to last about 10"* years. Even if 
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Right Ascension (J30l)0] 




Fig. 5 Upper part: images of various evolutionary phases of PWNe (from Gaensler & Slane [46]). 
From Left to right, up to down: X-ray image of the composite remnant G21. 5-0.9, free expansion 
phase; radio image of Vela SNR, displacement of the nebula due to the compression of the reverse 
shock during reverberation; SNR G 327.1-1.1 in radio (red) and X-rays (blue), relic PWN phase; 
W44 in radio, transition to the internal bow-shock phase; PSR B 1957+20 in Ha (green) an X-rays 
(red), ISM bow-shock nebula. Lower part: numerical hydrodynamical simulations of the various 
phases of the PWN-SNR evolution (from ||104.,23J ). Each figure of the lower part corresponds to 
systems shown in the upper one. 
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energy injection from the pulsar at these later times is negligible, PWNe can still 
be observed, due to the re-energization during compression. The interaction with 
the reverse shock can lead to a variety of different morphological structures if one 
consider also the pulsar proper motion Ill03lll04ll28ll43l . The most likely outcome 
of the interaction is that the nebula can be displaced with respect to the location of 
the pulsar At the beginning this might result into a system where the pulsar is not 
located at the center of the radio non thermal emission, analogous to what is ob- 
served in Vela. As the system evolves the reverse shock will completely displace the 
body of the PWN, creating a relic nebula. The relic PWN will mostly contain low 
energy particles, and will be visible in radio, while high energy particles, observ- 
able in X-rays will be seen only close to the pulsar G327.1 shows indeed this kind 
of morphology 1461 . Depending on projection effects one might also end up with 
very small X-ray nebulae centered on vast and large radio nebulae. In this regard 
polarization might prove essential to disentangle the structure. In particular, if the 
Rayleigh-Taylor instability is efficient and the radio nebula is disrupted and mixed 
with the ejecta, one expects to find a low level of polarization and no evidence for a 
global toroidal field. To this one must add the possibility that the ISM magnetic field 
might be dragged inside the PWN, or the nebular field inside the SNR shell. Given 
the importance of magnetic configurations for particle (cosmic ray) diffusion, one 
understands how important a proper model of older objects is. 

At later time the SNR ejecta starts cooling and the pulsar will eventually become 
supersonic. Once this happen the pulsar will form around itself a bow-shock PWN, 
and one expects an emission tail to form connecting the pulsar to the relic PWN. 
This model apply to the morphology and structure of W44 1461 . or possible IC443. 
Interestingly, the location with respect to the SNR where this happens does not 
depend on the pulsar proper motion, and turns out to be ^ 70% of the radius of 
the forward shock. An obvious question is if the jet-torus structure, that is observed 
in young and relatively undisturbed systems, can survive the later interaction with 
the SNR shell and the reverse shock. We know that in Vela the jet-torus is visible 
and does not appear to be distorted, this suggest that, as long as the evolution is 
subsonic, in the pulsar vicinity the dynamic of the nebula flow will still be regulated 
by the pulsar wind. On the other hand bow-shock simulations have shown that the 
typical size of the nebula in the head of the bow-shock and the flow dynamics do 
not allow the formation of collimated structures l24l l28l 11071 . In the case of SNR 
G327. 1-1.1 ll98l have suggested the presence of a jet-torus, but photons count and 
resolution are not high enough to make a definite statement. Deeper observation of 
transitional objects are needed. 

The ultimate phase of a PWN evolution depends on the pulsar kick. For slow 
moving pulsars the PWN will expands adiabatically inside the heated SNR, now in 
Sedov phase. Given the absence of energy injection, a PWN in this stage is probably 
only observable as a faint extended radio source, or possibly as a large TeV nebula 
due to IC from the relic leptons. To some extent this relic particle population might 
contribute to the diffuse gamma-ray background. On the contrary a fast moving 
pulsar can escape from the SNR, and will give rise to a bow-shock nebula due to 
the interaction with the ISM, through which it is moving at supersonic speeds 
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|28lll07| . These objects are observed both in Ha emission, due to ionization of ISM 
neutral hydrogen, and as long extended cometary-like source of non thermal radio 
and X-ray emission, due to the shocked pulsar wind, now forced to flow in the 
direction opposite to the pulsar motion. These nebulae might constitute one of the 
primary sources of positrons in the galaxy. 

Bow-shock PWNe, constitute a very interesting class among PWNe, and have re- 
cently received some attention, in particular regarding their X-ray emission. MHD 
models predict that the outflow in the tail of the bow-shock should have high speeds, 
of order of 0.5c, and that the tail should form a very well coUimated channel, with 
cross section comparable with the bow-shock size 1241 l28l [TOTl . However observa- 
tions have shown that in general the tail is wider than expected and that typical flow 
speed are high but of order 10000 km/s. It has been suggested that some form of 
mixing with the ISM, either via shear instability between the fast relativistic tail and 
the surrounding slower ISM, or via some particle contamination by ionized neutrals 
coming from the ISM, might be at play. A more detailed study of the fate of pairs 
injected in the tails of bow-shock nebulae, is essential, to assess the importance of 
pulsar as contributor to the pair CR background. 



4 Conclusion 

In the last few years, the combination of high resolution observations, and numerical 
simulations, has improved our understanding of the evolution and internal dynamics 
of PWNe. We can reproduce the observed jet-torus structure and we can relate the 
formation of the jet in the post shock flow to the wind magnetization. Simulated 
maps can reproduce many of the observed features, including the details of spectral 
properties. Results suggest that, the best agreement is achieved in the case of a 
wind with a large striped zone, even if MHD simulations are not able to distinguish 
between dissipation of the current sheet in the wind or at the TS. Results also suggest 
that it is possible to use X-rays imaging to constrain the pulsar wind properties; 
already the rings and tori observed in many PWNe have been used to determine the 
spin axis of the pulsar 1921 . Interestingly in Crab the inner ring appear less boosted 
in X-rays than the optical wisps (which shoud trace the same flow structure) are. 

Despite the undeniable successes of the MHD model, which are universally rec- 
ognized within the community, and the fact that there is general agreement that the 
observed X-ray properties are strongly dependent on the internal dynamic at the 
termination shock, for reasons unknown to the author, the old KC84 model is still 
used as a canonical reference for interpreting observations. While this might be un- 
derstandable for Radio or Optical data (where emission is quite homogeneous), it is 
completely unreasonable in X-ray. The fact that KC84 model is analytic and simple, 
is no excuse for its use, when it is clear that it is both qualitatively and quantitatively 
wrong, to the point that it basically fails to explain almost every single aspect of 
X-ray data. Moreover numerical tools and facilities are today widely available to 
conduct a correct study. 
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Interestingly, even when data are unreasonably averaged over spherical shells 
(even if imaging shows no hint of sphericity), they cannot be modeled using the 
original KC84 structure, and arbitrary velocity profiles are often assumed to repro- 
duce the data. So one trades a model (KC84) which is wrong (nebulae are far from 
spherical), but at least dynamically consistent (the correct solution of MHD equa- 
tions), for models which not only are wrong but also dynamically inconsistent, to 
the point that informations derived in this way have zero scientific valence, and the 
entire analisys is nothing more than a fit to the data. 

There are still however unsolved questions, and possible future developments for 
research in this field. All present simulations are axisymmetric, and none is able 
to address the problem of the stability of the toroidal field, nor can they reproduce 
the observed emission from the jet. It is not clear if small scale disordered field 
is present in the inner region: maybe a residual of the dissipation in the TS of the 
striped wind, or an outcome of the turbulence injected by the SASI-like instability 
of the TS. A combination of simulations and polarimetry might help to answer this 
question. 

Perhaps the more interesting developments might come from either the study of 
old systems, or from the investigation of particle energy distribution signatures in 
young ones. For the former, what is really needed is a large parameter study, where 
the various effects of pulsar proper motion and SNR reverberation are taken into 
account together with some simplified treatment of the spectral properties of these 
nebulae, in order to go beyond the simple qualitative morphological agreement, and 
provide templates for emission and spectral properties to be compared with obser- 
vations. For the latter, a study of possible signature of different injection mechanism 
at the TS, either as a function of shock properties, latitude or time, should be car- 
ried on by following the full particle distribution function in the nebula, instead of 
the simple power-law assumption. This will allow to verify if observable signature 
should be expected, in what band, and provide some constrain on the physics at the 
TS, and more important, on relativistic shock acceleration in general. 
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